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We report a measurement of the time-dependent CP asymmetry of —>■ D^^ph^ decays, where 
the light neutral hadron is a tt®, 77 or cu meson, and the neutral D meson is reconstructed in the 
CP eigenstates K^K~, Kgn^ or KgUJ. The measurement is performed combining the final data 
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samples collected at the T (dS") resonance by the BABAR and Belle experiments at the asymmetric- 
energy B factories PEP-II at SLAC and KEKB at KEK, respectively. The data samples contain 
(471±3) X 10® BB pairs recorded by the BABAR detector and (772±11) x 10® BB pairs recorded by the 
Belle detector. We measure the CP asymmetry parameters —rjfS = -|-0.66±0.10 (stat.)±0.06 (syst.) 
and C = —0.02 ± 0.07 (stat.) ± 0.03 (syst.). These results correspond to the first observation of CP 
violation in B® —>■ decays. The hypothesis of no mixing-induced CP violation is excluded in 

these decays at the level of 5.4 standard deviations. 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw 


In the standard model (SM) of electroweak interac¬ 
tions, CP violation arises from an irreducible complex 
phase in the three-family Cabibbo-Kobayashi-Maskawa 
(CKM) quark-mixing matrix [1]. The BaBAR and Belle 
experiments have established CP violating effects in the 
B meson system mg. Both experiments use their mea¬ 
surements of the mixing-induced CP violation in 6 —> ccs 
transitions m n] to determine precisely the parame¬ 
ter sin(2/3) = sin(2(/)i) [BaBAR uses (3 and Belle uses 
(j)i, hereinafter /3 is used]. The angle /3 is defined as 
arg [—where Vij is the CKM matrix ele¬ 
ment of quarks i, j. 

A complementary and theoretically clean approach to 
access /3 is provided by B^ —>• decays, where 

G {7r°,77,a;} denotes a light neutral hadron. These 
decays are dominated by CKM-favored b -A cud tree am¬ 
plitudes. CKM-disfavored b -A ucd amplitudes carrying 
different weak phases contribute also to the decays, but 
are suppressed by ~ *^-02 relative to the 

leading amplitudes. An interference between the decay 
amplitudes without and with B^-B^ mixing emerges if 
the neutral D meson decays to a CP eigenstate Dcp- Ne¬ 
glecting the suppressed amplitudes, the time evolution 
of B° -A D^ph° decays is governed by /3 [5] . Because 
only tree-level amplitudes contribute to B^ -A 
decays, these decays are not sensitive to most models of 
physics beyond the standard model (BSM). However, the 
measurement of the time-dependent CP violation enables 
testing of the measurements of 6 —> ccs transitions PIT] 
and provides a SM reference for the BSM searches in the 
mixing-induced CP violation of 6 —>■ s penguin-mediated 
B meson decays [SHU] ■ Any sizable deviation in the CP 
asymmetry of -A D^^phP decays from processes in¬ 
volving b -A CCS or penguin-mediated 5 —^ s transitions 
would point to BSM. Such deviations could for exam¬ 
ple be caused by unobserved heavy particles contributing 
to loop diagrams in 6 —?► ccs or 6 —?► s penguin transi¬ 
tions [T3] . 

An experimental difficulty in the use of B^ -A D^^phP 
decays arises from low B and D meson branching frac¬ 
tions [0(10“^) and 0{< 10“^), respectively] and low 
reconstruction efficiencies. Previous measurements per¬ 
formed separately by the BaBAR and Belle collaborations 
were not able to establish CP violation in these or related 
decays [TIHTB] . 

In this Letter, we present a measurement of the time- 


dependent CP violation in B^ -A D^^phP decays. For the 
first time, we combine the large final data samples col¬ 
lected by the BaBAR and Belle experiments. This new 
approach enables time-dependent CP violation measure¬ 
ments in the neutral B meson system with unprecedented 
sensitivity. 

The time-dependent rate of a neutral B meson decay¬ 
ing to a CP eigenstate is given by 
g-|At|/r^o 

g(At) = — - -{I -|- (7[iS sin(A?7idAt) 

Atbo 

-Ccos(AmdAf)]}, (1) 

where q = -|-1(—1) represents the 6-flavor content when 
the accompanying B meson is tagged as a B^ (i3°) and 
At denotes the proper time interval between the decays 
of the two B mesons produced in an T (45') decay. The 
neutral B meson lifetime is represented by Tgo, and the 
B^-B^ mixing frequency by Anid- Neglecting the CKM- 
disfavored decay amplitudes in —>■ D^^ph° decays, the 
SM predicts S = —r]f sin{2(3) and C = 0, where rjf is 
the CP eigenvalue of the final state, and S and C, re¬ 
spectively, quantify mixing-induced and direct CP viola¬ 
tion [Tt] . 

This analysis is based on data samples collected at 
the T(45) resonance containing (471 ± 3) x 10® Hi? 
pairs recorded with the BaBAR detector at the PEP-II 
asymmetric-energy e+e” (3.1 on 9 GeV) collider [TH] and 
(772±11) X 10® BB pairs recorded with the Belle detector 
at the KEKB asymmetric-energy e+e“ (3.5 on 8 GeV) 
collider [121. At BaBAR (Belle) the T{AS) is produced 
with a Lorentz boost of /?7 = 0.560 (0.425), allowing the 
measurement of At from the displacement of the decay 
vertices of the two B mesons. The BaBAR and Belle de¬ 
tectors are described in Refs. [2011^ . 

Reconstructed tracks of charged particles are consid¬ 
ered as kaon and pion candidates. Kaons are identi¬ 
fied using the particle identification techniques described 
in Refs. |20l [21]. Photons are reconstructed from en¬ 
ergy deposits in the electromagnetic calorimeters, and 
the energy of photon candidates is required to be at least 
30 MeV. Gombinations of two photons are considered as 
TT® meson candidates if the reconstructed invariant mass 
is between 115 and 150 MeV/c^. Gandidate rj mesons are 
reconstructed in the decay modes 77 —>■ 77 and 7r“'"7r“7r®. 
The invariant mass is required to be within 20 MeV/of 
the nominal mass [22] for 77 —> 77 candidates, and within 
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10 MeV/c^ for j] —)■ 7 r+ 7 r“ 7 r° candidates. For each photon 
in the ry —)■ 77 decay mode a minimal energy of 50 MeV 
is required. 

For uj mesons the decay mode lu —>■ 7 r“'" 7 r“ 7 r'^ is re¬ 
constructed with invariant mass required to be within 
15 MeV/c^ of the nominal mass [52]. Neutral kaons are 
reconstructed in the decay mode Kg —>■ tt+tt”, with in¬ 
variant mass required to be within 15 MeV/c^ of the 
nominal mass |22j . The requirements exploiting the Kg 
decay vertex displacement from the interaction point (IP) 
described in Refs. nulls] are applied. Neutral D mesons 
are reconstructed in the decay modes to CP eigenstates 
Dcp —> K'^K~, Kgir'^ and KgUj. The invariant mass is 
required to be within 12 MeV j(? of the nominal mass [22] 
for Dcp K^K~ candidates, within 25 MeV/c^ for 
Dcp —>■ KgTT^ candidates, and within 20 MeV/c^ for 
Dcp —>■ KgUj candidates. We reconstruct D*° mesons 
in the decay mode D*^ Zl° 7 r°, and the invariant mass 

must be within 3 MeV jc? of the nominal mass [22]. 

Neutral B mesons are reconstructed in the CP-even 
{rjf = -fl) final states P° —>■ Dcptt^ and Dcpr] (with 
Dcp —>■ Kgir^, KgUj), —)■ Dcpuj (with Dcp —>■ iFg 7 r°), 
B^ —>■ and DcpT] (with Dcp —>■ K^K~), and in 

the CP-odd (7/ = —1) final states P° —>■ Dcptt^, Dcprj, 
Dcpu> (with Dcp —>■ K^K~), and P° —>■ P^p 7 r° and 
D^pi] (with Dcp Kgir^) [51] . 

Neutral B mesons are selected by the beam- 
energy-constrained mass Mbc = toes = 
- (.Pb/c? [PaBAR uses Toes and Belle 
uses Mbc, hereinafter Mbc is used] and by the energy 
difference AE = Eg — Pbeanu where Pbeam denotes the 
energy of the beam, and pg and Eg are the momentum 
and energy of the B meson candidates, evaluated in 
the e“''e“ center-of-mass (c.m.) frame. The selected 
regions are 5.2 GeV/c^ < Mbc < 5.3 GeV/c^ and 
-100 MeV < AE < 100 MeV, except for B° Pcp7r° 
decays, where —75 MeV < AE < 100 MeV is required to 
exclude tails from partially reconstructed B~ —>■ D^*^^p~ 
decays peaking at AE « —250 MeV. 

In P° — D^uj and in ^ KgUj decays, the w vec¬ 
tor mesons are polarized. The angular distribution of 
u! —>■ 7 r“'" 7 r“ 7 r*^ decays is exploited to discriminate against 
background. The quantity cos Opr is defined as the co¬ 
sine of the angle between the neutral B meson direction 
and the normal to the 7 r"'' 7 r“ 7 r° plane in the w meson rest 
frame. A requirement of |cos0jv| > 0.3 is applied. 

After applying the above selection requirements, the 
average multiplicity of reconstructed P° —>■ D^^phP can¬ 
didates in an event is 1.3. In case of multiple B meson 
candidates in an event, one candidate is selected using 
a criterion based on the deviations of the reconstructed 
D^*'^ and ft.® meson masses from the nominal values. The 
probability for this method to select the correct signal is 
82% (81%) for BABm (Belle). 

In P® —>■ D^^ph^ decays, the dominant source of back¬ 


ground originates from e+e” —> <79 (g € {u,d, s,c}) con¬ 
tinuum events. This background is suppressed by using 
neural network (NN) multivariate classifiers that combine 
information characterizing the shape of an event [25j . 
The observables included in the NNs are the ratio R 2 
of the second to the zeroth order Fox-Wolfram moment, 
a combination of 16 modified Fox-Wolfram moments [IS] , 
the sphericity of the event [27] , and cos Og , where 9g is 
the angle between the direction of the reconstructed P 
meson and the beam direction in the c.m. frame. The 
NN selection reduces the background by 89.3% (91.8%) 
and has a signal efficiency of 75.5% (74.3%) for BaBAR 
(Belle). 



FIG. 1. The Mbc distributions (data points with error bars) 
and fit projections (solid line) of P® —>■ D^ph^ decays for (a) 
BABAR and (b) Belle. The dashed (dotted) lines represent 
projections of the signal (background) fit components. 

The signal yields are determined by unbinned maxi¬ 
mum likelihood fits to the Mbc distributions. In the fits, 
the signal component is parametrized by a Grystal Ball 
function |28j and the background component is modeled 
by an ARGUS function [29]. The experimental Mbc dis¬ 
tributions and fit projections are shown in Fig. [^ The 
signal yields are summarized in Table [^ 


TABLE I. Summary of P® —^ P^ft® signal yields. 


Decay mode 

BABAR 

Belle 

P® ^ PcpTT® 

241 ± 22 

345 ± 25 

P® —>■ DqpT] 

106 ± 14 

148 ± 18 

P® — >• Dcpijj 

66± 10 

151 ± 17 

P® PcpTT® 

72± 12 

80± 14 

P® — >■ DcpP 

39 ±8 

39± 10 

P® ^ P^^ft® total 

508 ± 31 

757 ± 44 


The time-dependent CP violation measurement is per¬ 
formed using established BABAR and Belle techniques for 
the vertex reconstruction, the flavor-tagging, and the 
modeling of At resolution effects (see Refs. [51171 1501155] 1 
and is briefly summarized below. The proper time inter¬ 
val At is given as where Az is the distance between 
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the decay vertices of the signal B meson and of the ac¬ 
companying B meson. The B^ —>■ D^^phP signal decay 
vertex is reconstructed by a kinematic fit including in¬ 
formation about the IP position. For Belle, an iterative 
hierarchical vertex reconstruction algorithm following a 
bottom-up approach starting with the final state particles 
is applied, while for BaBAR the vertex reconstruction in¬ 
cludes simultaneously the complete B meson decay tree 
including all secondary decays. In the kinematic fits, 
the invariant masses of 7r°, ij, uj, and Dcp candidates 
are constrained to their nominal values |22j . The decay 
vertex and the 6-flavor content of the accompanying B 
meson are estimated from reconstructed decay products 
not assigned to the signal B meson. The 6-flavor con¬ 
tent is inferred by flavor-tagging procedures described in 
Refs. [5115^. The applied algorithms account for different 
signatures such as the presence and properties of prompt 
leptons, charged kaons and pions originating from the 
decay of the accompanying B meson, and assign a flavor 
and an associated probability. Selection requirements on 
the quality of the reconstructed decay vertices and the 
At measurements are applied. 

The CP violation measurement is performed by maxi¬ 
mizing the log-likelihood function 

In £ = ^ In ’ (2) 

i 3 

where the indices i and j denote the events reconstructed 
from BaBAR and Belle data, respectively. The probability 
density function (p.d.f.) describing the At distribution 
for BaBAR is defined by 

k 

and for Belle by 

pBelle ^ 5] /fe / [Pk {At') Rk {At - At')] d {At') 

k 

+/ 01 P 0 I (At), (4) 

where the index k represents the signal and background 
p.d.f. components. The symbol P^ denotes the p.d.f. de¬ 
scribing the proper time interval of the particular phys¬ 
ical process, and Rk refers to the corresponding reso¬ 
lution function. The fractions fk are evaluated on an 
event-by-event basis as a function of Mbc. Belle treats 
outlier events with large At using a broad Gaussian func¬ 
tion in the p.d.f. component PqI with a small fraction of 
/oi « 2 X 10“"^, while BaBAR includes outlier effects in the 
resolution function. The signal p.d.f. is constructed from 
the decay rate in Eq. Q , including the effect of incorrect 
flavor assignments and convolution with resolution func¬ 
tions to account for the finite vertex resolution. The mod¬ 
els of the At resolution effects at BaBAR and Belle follow 
different empirical approaches and are described in detail 


in Refs. The background p.d.f.s for BaBAR and 

Belle are composed of the sum of a Dirac delta function 
to model prompt background decays and an exponential 
p.d.f. for decays with effective lifetimes. The background 
p.d.f. is convolved with a resolution function modeled as 
the sum of two Gaussian functions. The background pa¬ 
rameters are fixed to values obtained by fits to the events 
in the Mbc < 5.26 GeV/c^ sidebands. 



FIG. 2. (color online). The proper time interval distributions 
(data points with error bars) for R® tags (red) and tags 
(blue) and the CP asymmetries of B° —> D^^phP decays for 
(a)-(b) BABAR and (c)-(d) Belle for candidates associated with 
high quality flavor tags. The solid lines show projections of 
the sum of signal and background components in the fit, while 
the hatched areas show only the background components. 

The combined BaBAR and Belle measurement is 
performed by maximizing Eq. (1^ for events in the 
5.27 GeV/c^ < Mbc < 5.29 GeV^ signal region. The 
values of rgo and Auid are fixed to the world aver¬ 
ages [22]. The free parameters in the fit are S and C. 
The result is 

—i]fS = -1-0.66 ± 0.10 (stat.) ± 0.06 (syst.), 

C = -0.02 ± 0.07 (stat.) ± 0.03 (syst.). (5) 

The linear correlation between —T]fS and C is —4.9%. 
Through comparison of the log-likelihood of the fit to the 
































distribution from an ensemble test performed with input 
from the data distributions, a p-value of 0.46 is obtained. 
The flavor-tagged proper time interval distributions and 
projections of the fit are shown in Fig. 

The evaluation of the systematic uncertainties in the 
CP violation parameters follows standard approaches of 
the BaBAR and Belle experiments described in detail in 
Refs. 0 13 |33]; the results are summarized in Table [TI] 
For the vertex reconstruction, the sources of systematic 
uncertainties include the applied constraints and selec¬ 
tion requirements on the vertex fits of the signal B me¬ 
son and the accompanying B meson, and on the At fit 
range. These contributions are estimated by variations 
of the constraints and selection requirements. The sys¬ 
tematic uncertainties due to the misalignment of the sili¬ 
con vertex detectors are estimated by Monte Carlo (MC) 
simulations. For BAlBAR, the uncertainty of the z scale 
is estimated by variations of the z scale and correspond¬ 
ing uncertainties. For Belle, a possible At bias is esti¬ 
mated using MC simulations. The systematic uncertain¬ 
ties due to the At resolution functions, the parameter¬ 
ization of the At background p.d.f., the calculation of 
the signal purity, the flavor-tagging, and the physics pa¬ 
rameters tbo and Am^ are estimated by variation of the 
fixed parameters within their uncertainties. Fit biases 
are estimated using large samples of MC-simulated sig¬ 
nal decays. The contribution of backgrounds that have 
the same final states as the reconstructed B^ -A 
decay modes and that can peak in the Mbc signal region 
is estimated using D meson mass sidebands on data and 
using generic BB MC samples. These backgrounds ac¬ 
count for less than 8 % of the signal and consist mainly 
of flavor-specific decays such as partially reconstructed 
B~ -A decays. The systematic uncertainty due 

to this peaking background is estimated using MC simu¬ 
lations in which the peaking background is modeled, and 
the nominal fit procedure, which neglects this peaking 
background, is applied. The effect of interference be¬ 
tween b -A cud and b -A ucd decay amplitudes of the 
accompanying B meson is estimated using MC simula¬ 
tions that account for possible deviations from the time 
evolution described by Eq. 0 m- Possible correlations 
between BaBAR and Belle are accounted for in the evalu¬ 
ation of the contributions due to the physics parameters, 
the peaking background, and the tag-side interference. 
In the MC studies described above, the largest devia¬ 
tions are assigned as systematic uncertainties. The total 
systematic uncertainty is the quadratic sum of all contri¬ 
butions. 

The statistical significance of the results is estimated 
using a likelihood-ratio approach by computing the 
change in 2 In £ when the CP violation parameters are 
fixed to zero. The effect of systematic uncertainties is 
included by convolution of the likelihood distributions. 
No significant direct CP violation is observed. The mea¬ 
surement excludes the hypothesis of no mixing-induced 


TABLE II. Summary of systematic uncertainties for the time- 
dependent CP violation measurement in B° -A O'^pP’ decays 
(in units of 10“^). 


Source 

5 

C 

Vertex reconstruction 

1.5 

1.4 

At resolution functions 

2.0 

0.4 

Background At PDFs 

0.4 

0.1 

Signal purity 

0.6 

0.3 

Flavor-tagging 

0.3 

0.3 

Physics parameters 

0.2 

< 0.1 

Possible fit bias 

0.6 

0.8 

Peaking background 

4.9 

0.9 

Tag-side interference 

0.1 

1.4 

Total 

5.6 

2.5 


CP violation in B^ -A D^phP decays at a confidence level 
of 1 — 6.6 X 10“®, corresponding to a significance of 5.4 
standard deviations. 

The analysis is validated by a variety of cross-checks. 
The same measurement is performed for B'^ -A 
decays with the CKM-favored -A K~tt^ decay mode. 

These decays provide a kinematically similar, high statis¬ 
tics control sample. The result agrees with the assump¬ 
tion of negligible CP violation for these decays. Mea¬ 
surements of the neutral B meson lifetime using the 
control sample and B° -A D^pP decays yield Tgo = 
1.518±0.026 (stat.) ps and Tpo = 1.520±0.064 (stat.) ps, 
respectively, in agreement with the world average Tgo = 
1.519 ± 0.005 ps [ 33 . All measurements for the control 
sample and for B^ -A D^pP decays have also been per¬ 
formed for data separated by experiment and by decay 
mode, and yield consistent results. The results for B^ -A 
D^pP decays separated by experiment are sin(2/3) = 
0.52 ± 0.15 (stat.) for BaBAR and 0.83 ± 0.15 (stat.) for 
Belle, and the results separated by the CP content of the 
final states are sin(2/3) = 0.52 ± 0.15 (stat.) for CP-even 
and 0.80 ± 0.15 (stat.) for CP-odd. 

In summary, we combine the final BaBAR and Belle 
data samples, totaling more than 1 ab~^ collected at the 
T (4S) resonance [13 I3S], and perform a simultaneous 
analysis of the data collected by both experiments. We 
observe for the first time CP violation in P° -A D^pP 
decays driven by mixing-induced CP violation. We mea¬ 
sure sin(2/3) = 0.66±0.10 (stat.)±0.06 (syst.). This result 
agrees within 0.2 standard deviations with the world av¬ 
erage of sin(2/3) = 0.68±0.02 [53 measured from b -A ccs 
transitions, and is consistent with the measurements of 
6 —>■ s penguin-mediated B meson decays mn] at cur¬ 
rent precision. The presented measurement supersedes 
the previous BaBAR result for P° —>■ D^^pP decays [15] . 
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